Precambrian Research 409 (2024) 107453 


J 


ELSEVIER 


Contents lists available at ScienceDirect 
Precambrian Research 


journal homepage: www.elsevier.com/locate/precamres 


Check for 


Hydrothermal seawater eutrophication triggered local macrobiological upaaies 
experimentation in the 2100 Ma Paleoproterozoic Francevillian sub-basin 


Ernest Chi Fru’, Jérémie Aubineau’‘, Olabode Bankole ©, Mohamed Ghnahalla €, 


Landry Soh Tamehe“, Abderrazak El Albani‘ 


* College of Physical and Engineering Sciences, School of Earth and Ocean Sciences, Centre for Geobiology and Geochemistry, Cardiff University, Cardiff CF10 3AT, 


Wales, UK 


> Géosciences Environnement Toulouse, CNRS UMR 5563 (CNRS/UPS/IRD/CNES), Université de Toulouse, Observatoire Midi-Pyrénées, Toulouse, France 
€ Université de Poitiers UMR 7285-CNRS, Institut de Chimie des Milieux et Matériaux de Poitiers, 5, rue Albert Turpin (Bat B35), 86073 Poitiers Cedex, France 


4 China Nonferrous Metals (Guilin) Geology and Mining Co. Ltd, Guilin 541004, China 


ARTICLE INFO ABSTRACT 


Keywords: 

Manganese carbonates 
Iron carbonates 
Proterozoic 

Congo Craton 
Seawater oxygenation 


It is thought that the global predominance of small-size unicellular prokaryotic life in the oceans until the 
emergence of large-size multicellular organisms to ecological dominance in the Ediacaran Ocean after 635 
million years ago (Ma), was partly constrained by paleo-dynamic nutrient limitation, with phosphorus (P) being 
the principal limiting resource. Here we couple an episode of intense submarine hydrothermal alteration of a 
nutrient-rich seafloor reservoir to the collision of the Congo-Sao Francisco cratons at ~2100 Ma, to unravel a 


paleo-geodynamic incident of seawater P enrichment in the Paleoproterozoic Francevillian sub-basin. We pro- 
pose that this previously unrecognized local pulse in dissolved seawater P concentration, of comparable 
magnitude to Ediacaran seawater levels, set the stage for Earth’s earliest biospheric experimentation towards 
macrobiological complexity ~2100 million years ago. 


1. Introduction 


Being a critical biolimiting nutrient over geological timescales 
(Tyrrell, 1998), low seawater phosphorus (P) concentration is frequently 
advanced as a key cause of mildly productive Proterozoic oceans, 
reduced organic carbon (Corg) burial rates and limited ocean—atmo- 
sphere oxygenation (e.g., Reinhard et al., 2017; Guilbaud et al., 2020; 
Laakso et al., 2020; Poulton et al., 2021; Alcott et al., 2022; Chi Fru 
et al., 2023; Dodd et al., 2023). This view is anchored on the significance 
of P being the backbone of nucleic acids, cell membrane phospholipids 
and cell energy trafficking by adenosine triphosphate (ATP). 

Thus, environmental phosphate (POs) bioavailability is thought to 
have exerted considerable influence on the pace of biological evolution 
from small unicellular life forms to large multicellular metazoans (e.g., 
Brocks et al., 2017; Reinhard et al., 2017, 2020; Eckford-Soper et al., 
2022). Further, the dependence of life on oxidoreductase enzymes 
containing reactive metal cores to generate cellular energy (Jelen et al., 
2016; Dupont et al., 2006, 2010; Moore et al., 2017), has spawned the 
hypothesis that fluctuations in trace metal abundance through Earth 
history masterminded their selective use in enzymatic reactions through 
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time (Saito et al., 2003; Dupont et al., 2006, 2010; Robbins et al., 2013; 
Scott et al., 2013; Moore et al., 2017). 

According to conventional thinking, unequivocal evidence for 
eukaryotic fossils first appeared in the geological record some 
1700-1600 million years ago (Ma), almost a billion years before the 
appearance of animals after 635 Ma (Knoll et al., 2006; Javaux, 2011; 
Buffield, 2009, Butterfield, 2015; Zhu et al., 2016; Reinhard et al., 2020; 
Eckford-Soper et al., 2022). Even more profound, some four hundred 
millions earlier, the ~2100 Ma Paleoproterozoic Francevillian sub-basin 
in southeastern Gabon (Fig. 1A), records a unique assemblage of enig- 
matic macroscopic specimens interpreted by some studies (e.g., El 
Albani et al., 2010, 2014, 2019, 2023; Ossa Ossa et al., 2023) and 
disputed by others (e.g., Anderson et al., 2016; Fakhraee et al., 2023), to 
represent models of Earth’s earliest oxygen-respiring colonial life forms, 
some which were capable of complex motility behaviour. Moreover, 
persistent association of the Francevillian macrofossils with a rich 
benthic microbial community (Aubineau et al., 2018, 2019, 2021), is 
consistent with the requirements for a high energy resource-rich 
ecosystem to enable evolutionary gigantism and complexity (Ward 
et al., 2012). Thus, the recent proposition of a macro-eukaryotic lifestyle 
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resulting in cellular Zn, Ni and Co enrichment in the Francevillian 
macrofossils (El Albani et al., 2023; Ossa Ossa et al., 2023) raises the 
question of whether a shift in environmental nutrient supply to the 
paleo-Francevillian sub-basin drove the emergence of these enigmatic 
macrofossils. 

In the modern ocean, small phytoplankton and picoeukaryotes 
dominate low nutrient habitats and larger complex organisms more 
productive nutrient-rich ecosystems. Similarly, it is thought that 
ecological dominance of higher eukaryotic algae in post-Cryogenian 
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Snowball PO}-rich seawater powered complex food webs that sus- 
tained more efficient nutrient and energy exchange to drive oxygen- 
dependent Ediacaran ecosystems towards larger and more complex 
metazoans (Brocks et al., 2017; Van Maldegem et al., 2019). Thus, 
environmental oxygen concentration increasing to above a certain 
minimum threshold may likely have been a key contributor to greater 
macrobiological complexity (e.g., Planavsky et al., 2014; Cole et al., 
2020; Reinhard et al., 2020; Eckford-Soper et al., 2022; Craig et al., 
2023). Although the link between environmental oxygen threshold 
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Fig. 1. A) Geological map of the Francevillian Group showing location of the four sampled drill holes (KO42039, KO54040, KO42028 and KO65035). The Bangombe 
Plateau quarry that records the Francevillian biota is indicated by Mou. B) Lithostratigraphy of individual drill cores. C) Combined generalized drill core stratigraphy. 
FA to FE represents Franceville A, B, C, D and E Formations. The FB Formation is divided into FB1 and FB2 members composed of FB1a, FB1p and FB,, and FB, and 
FB units, respectively. The horizontal purple line is the reference datum for horizontalization, as represented by Reynaud et al. (2018), which allowed stratigraphic 


correlations between drill cores. 
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levels and the role of genetic innovation in the transition to complex 
multicellular life is debated (e.g., Buffield, 2009, Erwin et al., 2011; 
Mills et al., 2014; Bozdag et al., 2021), oxygen availability and genetic 
innovation, in tandem with nutrient levels, likely worked in synergy to 
drive biological diversification towards greater complexity, given that 
low environmental oxygen levels also select for small size individuals 
and higher levels large size organisms (Sperling et al., 2013). 

In line with this view, there is compelling evidence linking a stable 
shallow oxygenated continental seawater margin setting with emer- 
gence of the Francevillian macrofossils (El Albani et al., 2010, 2014). 
This is intriguing because the Francevillian macrofossils appear on the 
heels of the 2450-2200 Ma Great Oxidation Event (GOE) (Poulton et al., 
2021) as are the Ediacaran macrofossils that coincide with the Neo- 
proterozoic Oxygenation Event (NOE). However, the role of nutrient 
bioavailability, and particularly POF, remains to be associated with 
gigantism in the Francevillian macrofossils. This has resulted in 
considerable uncertainties on what triggered their emergence and 
whether their biogenicity can be pinned to environmental conditions 
that enable macrobiological evolution. 

Here we explore the hypothesis that dramatic environmental change 
resulting in high nutrient bioavailability in seawater drives trophic level 
complexity and evolutionary gigantism (Ward et al., 2012; Reinhard 
et al., 2016), to explore the biogenic bridge between environmental 
change, seawater nutrient bioavailability and the appearance of mac- 
roorganisms in the Francevillian sub-basin at ~2100 Ma. To achieve this 
goal, we undertake a systematic reconstruction of a whole suite of 
seawater biogeochemical conditions that prevailed just before the 
Francevillian macrofossils that grew to >17 cm long forms, some as 
lenticular, segmented, tubular, lobate and string shaped assemblages (El 
Albani et al., 2010, 2014, 2019, 2023; Ikouanga et al., 2023), emerged. 

The biogeochemical evidence surprisingly reveals unprecedented 
transitory enrichment of the ~2100 Ma Paleoproterozoic Francevillian 
sub-basin seawater margin with P levels comparable to those that sup- 
ported the rise of the Ediacaran biota just before the appearance of the 
Francevillian macrofossils. We track the trigger for this abrupt local 
environmental change to hydrothermal alteration of a nutrient-rich 
seafloor reservoir initiated by collision of the Congo-Sao Francisco cra- 
tons during the Eburnean-Transamazonian orogeny. 


2. Geological setting 


The up to 2.5 km-thick epicontinental sedimentary package of the 
~2100 to 2000 Ma Francevillian Group in southeastern Gabon and its 
exquisitely preserved and unmetamorphosed rocks, comprise five 
distinct formations emplaced unconformably on an Archean granitoid 
and greenstone belt basement (Fig. 1A). The ~1000 m-thick basal FA 
Formation, mainly of fluvial to deltaic sandstones, is rich in organic 
matter, hematite and uranium. 

The overlying FB, FC, FD, FE formations were deposited in pre- 
dominantly marine environments (Mossmann, et al., 2005; Canfield 
et al., 2013; Reynaud et al., 2018), with the FB,.. and FB2,4 units 
dominated by fine-grained shales containing a Mn-rich ore (protore) 
made up of Mn-rich carbonates underlain by a thin Mn-rich sandstone 
layer and an Fe-rich formation (Fig. 1A) composed mainly of Fe-rich 
carbonates (Canfield et al., 2013; Ngombi Pemba et al., 2014; Rey- 
naud et al., 2018). The FB1b unit comprising heterolithic beds (banded 
shales), together with the Fe-Mn-rich lithologies of the FBlc unit, 
constitute the lithologies analysed in this study (Fig. 1A). 

The FB2 Member hosts the macrofossils, which are associated with 
numerous microbial mat-related structures (Aubineau et al., 2018; El 
Albani et al., 2010, 2014, 2019; Reynaud et al., 2018), while Mn 
deposition in the upper FB; Member is believed to correlate temporally 
with deposits in the Sao Francisco Basin, Brazil, at the time amalgam- 
ated to the Francevillian Group (dos Santos et al., 2023). Correlation of 
8!3Cearb and mene in the FB and lower FC formations with the Zaonega 
Formation, Russia, together with geochronology, constrain the 
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depositional age of the Francevillian Group to ~2.11-2.06 Ga. (Bros 
et al., 1992; Thiéblemont et al., 2009; Kump et al., 2011; Canfield et al., 
2013; Sawaki et al., 2017). 


3. Materials and methods 
3.1. Sample collection and preparation 


Stratigraphic reconstruction relied on identification of the Fe-rich 
carbonate facies as a transitional boundary marker from the FBıp to 
the FB;, unit (Weber, 1968; Gauthier-Lafaye and Weber, 2003; Canfield 
et al., 2013; Reynaud et al., 2018). The latter studies have extensively 
described the sedimentological structures and sequence stratigraphic 
framework of the FB Formation. This well-recognized Fe-rich carbonate 
facies described in Canfield et al. (2013) enabled correlation between 
drillholes (KO42028, KO65035, KO54040 and KO42039) obtained from 
the Okouma Plateau area (Fig. 1B, C). The drillholes intersect various 
sections, thickness, and depths of the FBjp.. units (Reynaud et al., 2018), 
with the deepest KO42028 reaching ~325 m and KO65035, KO54040 
and KO42039 ~150, 140 and 150 m, respectively. Four main lithologies 
common to the drillholes include the 51-240 m thick laminated heter- 
olithic siltstones (banded shales), characterized by creamy-whitish 
dolomitic bands alternating with siltstones and greyish shales 
(Figs. 1C, D & 2A-F). These are overlain by 5-10 m thick pyritized 
shales, and then by an 8-17 m thick Fe-rich carbonate facies, sometimes 
containing breccias and pyrite laminae. This passes over to a thin pyri- 
tized shale unit, overlain by 41-53 m thick Mn-rich black shales with a 
thin Mn-rich sandstone horizon at the base. 


3.2. Petrography and mineralogy 


Polished thin sections were prepared for representative samples and 
examined under a polarized Nikon Eclipse E-600 optical microscope at 
the University of Poitiers, France. The thin sections were carbon-coated 
and studied with a JEOL scanning electron microscope (SEM) coupled 
with energy dispersive X-ray (EDX) for textural, mineralogical and spot 
elemental composition in backscattered electron (BSE) mode. X-ray 
diffraction (XRD) on bulk rock and <2 um clay fractions was conducted 
with a Bruker D8 ADVANCE diffractometer as previously established by 
Moore and Reynolds (1997) and Ngombi Pemba et al. (2014). 


3.3. Bulk rock geochemistry 


Approximately 300 mg of powdered samples were fused with 900 mg 
ultra-pure lithium metaborate (LiBO2) for trace and major element an- 
alyses (Carignan et al., 2001; Ngombi Pemba et al., 2014). Major ele- 
ments (MEs) were measured with an inductively coupled plasma optical 
emission spectrometer (ICP-OES) and trace elements (TEs) by ICP mass 
spectrometry (ICP-MS) at Service d'Analyse des Roches et Minéraux 
(SARM), Nancy, France. Measurement uncertainties and detection limits 
are presented in Tables S1 and S2. Analytical procedures were validated 
with basalt BR, anorthosite AN-G, diorite DR-N, serpentinite UB-N and 
granite GH as geological reference standards. Selected redox sensitive 
metals, major elements and TEs were normalized with Al as a detrital 
indicator and standardized to Post-Archean Australian Shale (PAAS) 
composition (Taylor and McLennan, 1985) to distinguish authigenic 
from detrital enrichments. Elemental enrichment factor (EF) was 
calculated as EFx = [(X/Al sample/(X/ADpaas], where X = concentration 
of element of interest (Tribovillard et al., 2006). 


3.4. Chemical index of alteration (CIA) 


The chemical index of alteration (CIA) was used to quantify the de- 
gree of chemical weathering of source rocks, where CIA = 
[Al203/(Al203 +CaO* +Na20 +K20)] x 100, in molar proportions (Fedo 
and Babechuk, 2023). For the CIA analysis, CaO* represents the quantity 


E. Chi Fru et al. 


G 


Fe-dolomite 


Pyrite 
Clays 


x 
k 


Pyrite 


J 


Dolomite ° a" 
| Quartz 


rhs Ankerite 


Pyrite 


Precambrian Research 409 (2024) 107453 


Mnerich black shale, 36.8 m 
Mp-rich sandstone, 69.4 m 
Fe-rich carbonate, 75.8 m 
Pyritized black shale, 72.9 m 


Fe-rich carbonate, 77.4 m 


Pyritized black shales, 80.95 m 


Banded shale, 117.6 m 


Banded shale, 229.3 m 


Banded shale, 322.25 m 


— a 
— EG 
Mn-rich black shale, 36.8 m 


r 


f 


Mn-tich black shale, 44.7 m 


Mn-tich sandstone, 69.4 m 


Pyritized black shale, 73.25 m 


Fe-rich carbonate, 77.4 m 


Pyritized black shale, 80.95 m 


Banded shale, 162.5 m 


oH 
Y 


Banded shale, 255.7 m 


Banded shale, 279.1 m 


Banded shale, 318.6 m 


Banded shale, 322.25 m 


Fig. 2. Lithofacies, petrographic and mineralogical analysis of representative lithologies for the studied section. A) Mn-rich black shale. B) Mn-rich sandstone. C) 
Pyritized black shale. D) Fe-rich carbonate shales. E) Dolomitic and ankeritic facies showing breccias composed of dolomitic clasts. F) Banded shales showing white 
bands of silty dolomite alternating with black shale units. G-H) Dispersed pyrite engulfing clays and Fe-dolomite within the Fe-rich carbonate. I-J) Dolomitic and 
ankeritic minerals in the Mn-rich black shales. K-L) Clay and ankerite minerals filling the matrix of the underlying banded shale. M) Bulk rock X-Ray Diffractograms 
(XRD) for representative samples. N) XRD spectra for <2 um clay fraction for representative samples from the Okouma Plateau. Black spectra are for air-dry (AD) 
samples and red for ethylene-glycol (EG). Chl, chlorite; I/M, illite/micas; I/S, illite/smectite; Alb, albite; Q, quartz; Py, pyrite; An, ankerite, Do, dolomite; Rh, rhodochrosite. 
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of CaO associated with silicate phases after correction for CaO-rich 
phases (McLennan et al., 1993). We assumed moles of silicate CaO* = 
Na2O when CaO > Na2O in molar proportions and CaO* = CaO when 
CaO < Na2O in molar proportions (e.g., Jian et al., 2013). Weathering 
changes and sediment source were further inferred by projecting the 
obtained CIA values on an Al203-(CaO* + Na20)-K20 ternary diagram. 
CIA values >80 depict intense leaching of mobile Ca**, Na‘, K* cations 
relative to immobile Al** during chemical weathering of parent rock, 
while values <50 represent weak to near absence of chemical 
weathering. 

The CIA values of sedimentary rocks can be influenced by input 
source composition, chemical weathering, sediment recycling, sorting 
and diagenesis, which can complicate comparisons of sedimentary ar- 
chives (Fedo et al., 1995; Fedo and Babechuk, 2023). We therefore 
considered the effects of these processes on our CIA values. The data 
were interpreted as CIA < 50 % for unweathered igneous rocks, 60-80 % 
for moderately weathered rocks and >80 % for highly weathered rocks 
(Fedo et al., 1995; Fedo and Babechuk, 2023). 


3.5. Redox reconstruction 


Rare Earth Elements (REEs), redox sensitive metals, C and N sys- 
tematics and Fe speciation were combined to assess redox depositional 
conditions on 58 representative samples, bracketing the FB1p-c units. 
Iron speciation was performed at Cardiff University, UK, following the 
refined Fe extraction method of Poulton and Canfield (2005), taking into 


consideration pitfalls caused by elevated carbonate and low Fe content, 
as reviewed by Raiswell et al. (2018). Data were compared with previ- 
ously published results for the correlative Bangombe Plateau (Canfield 
et al., 2013). Instrument interference on Eu composition caused by Ba 
was evaluated as described in Fig. S1 and shown to be insignificant. 


3.6. Carbon and nitrogen stable isotope analysis 


Organic carbon (6'3Corg) and whole rock 5'°Npuk isotope analyses 
were performed on the same samples used for Fe-based redox recon- 
struction. First, powdered samples were treated with 6 N HCl at 70 °C for 
1 h to remove carbonate minerals. After repeated washing with distilled 
water and centrifugation to neutral pH, the residue was dried in a clean 
hood overnight. Small aliquots of de-carbonated and whole-rock 
powdered samples were weighed and sealed in tin capsules for C and 
N isotope analyses, respectively, with a blank inserted between mea- 
surements. Data were obtained by flash combustion on a Thermo Sci- 
entific DELTA V Advantage isotope ratio mass spectrometer (IRMS), 
operated under continuous helium flow mode, and coupled to a COST- 
ECH 4010 elemental analyzer, at the University of California, Riverside, 
USA. The data are reported in delta (6) notations, as part per mil (%bo) 
deviation relative to Vienna Pee Dee belemnite (VPDB) for carbon and 
atmospheric N3 for nitrogen. The 5'°Corg and 6!°Npuik data were cali- 
brated against two internal laboratory standards — acetanilide from two 
different batches (6'°C = —27.9 %o, 5!°N = —0.8 %o and 5!°C = —33.7 %o, 
5'°N = —0.8 %o) and Hawaiian glycine (6'°C = —36.6 %o, 5!°N = 11.3 %o) 
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and two international reference standards, SDO-1 (813°C = —30.0 %o, 
5'°N = —0.8 %o) and SGR (6!'°C = —29.3 %o, 5!°N = —17.4 %o) (Dennen 
et al., 2006). Replicate analyses produced a 0.09 %o standard deviation 
(10) for 5 Corg (n = 9) and 0.11 %o for 5 °Nbulk (n = 8). Total organic 
carbon (TOC) and total nitrogen (TN) concentrations were calculated 
from peak areas of CO2 and Ng, respectively. Reported 5 >Nbulk values 
are only for samples with TN concentrations above instrument detection 
limit. 


4. Results 
4.1. Lithostratigraphic enrichment of major and trace elements 


Consistent with lithostratigraphy (Fig. 2A-F) and mineralogical an- 
alyses (Fig. 2G-N), significant enrichment of a range of MEs, TEs and 
REEs, delineate a sharp transition from the underlying banded shales to 
the overlying Fe-Mn-rich lithologies (Fig. 3A-P; Tables S3-5). Particu- 
larly, considerable P205 (P) and Fe203 (Fe) enrichment differentiates the 
Fe-rich carbonates from the rest of the lithologies (Fig. 3A, M, N). The 
banded shales are enriched in detrital quartz, albite and chlorite, with 
illite/mica present to a lesser extent, while the Fe-Mn-rich lithologies 
containing siderite (FeCO3) and rhodochrosite (MnCO3) as dominant Fe 
and Mn minerals, respectively (Fig. 2M, N), exhibit negligible detrital 
incorporation, as are the pyritized black shales. The latter observations 
are consistent with normalization of TE content with detrital Al and with 
up sequence switch towards greater chemical weathering CIA values 
(Figs. 4, 5A & Table S3). 

The Al2O3 and TiOz content show similar co-enrichment trends 
across the section (Fig. 5B, C). Fe2O3 peaks in the Fe-rich carbonates and 
is within a comparable range in the Mn-rich black shales and banded 
shales (Fig. 5D). MnO concentrations show a systematic increase up 
section (Fig. 5E), with corresponding Fe203/Mn0O ratios decreasing in a 
similar direction (Fig. 5F). The ratios of Fez03/TiO2 and Fe203/Al,03 
(Fig. 5G, H) display a comparable behaviour as are those of MnO/TiO2 
and MnO/A1,03 ratios (Fig. 5I, J). Considering Al and Ti as a detrital 
tracer, the Fe203 content in the Mn-rich black shales and banded shales 
appears to have experienced greater dilution by detrital input compared 
to Fe203 in the Fe-rich carbonates (Fig. 5G, H). Similarly, detrital input 
is suggested to have exerted a stronger influence on MnO enrichment in 
the Mn-rich black shales and banded shales than in the Fe-rich carbon- 
ates (Fig. 5I, J). 

Except for a moderate positive correlation of 54 % for the banded 
shales, a weak linear correlation of 18 % is observed when Fe203/TiO2 
ratios for the Mn-rich black shales and Fe-rich carbonates are plotted 
against Fe2O3 concentrations (Fig. 6A)—being up to 71 % linearly 
correlated cross all three lithologies combined (Fig. 6B). On the other 
hand, MnO/TiO, ratios for the banded shales, Fe-rich carbonates and the 
Mn-rich black shales lithologies show a strong positive correlation of 96 
%, 82 % and 74 %, with MnO concentrations, respectively (Fig. 6C)— 
being only 47 % linearly correlated for all three compared (Fig. 6D). 


4.2. Indicators of hydrothermal activity 


The corresponding Fe/Ti versus Al/(Al + Fe + Mn) values for the 
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individual facies indicate that the Fe-rich carbonates were deposited 
under more influential hydrothermal activity, with one sample plotting 
in the intense hydrothermal field (Fig. 7A). While both the Mn-rich 
carbonates and the banded shales show lower Fe/Ti ratios, the latter 
lithologies are distinguished by extremely high Al/(Al + Fe + Mn) 
values. 

Overall, the samples group with banded iron formations (BIFs) of 
comparable age and of low temperature hydrothermal origin in the 
northwestern Congo Craton (Fig. 7B, C). However, cross plots of their 
(Sm/Yb)cyn versus (Eu/Sm)cy ratios (Fig. 7C) indicate subtle geochem- 
ical differences consistent with varying degrees of hydrothermal fluid 
composition and seawater mixtures. Distinctively, four of the Fe-rich 
carbonates and one each for the pyritized and Mn-rich black shales 
fall between 0.1-1 % hydrothermal fluid composition and up to 1-5 % 
for one Fe-rich carbonate sample (Fig. 7C). Although less evident on the 
Y/Ho versus Eu/Eu* plot (Fig. S2A), hydrothermal fluid contribution to 
the facies is further indicated by PAAS normalized REE spider diagrams 
with variable positive Eu anomalies (Fig. 7D-G). 

Moving up sequence, the sediments record average Y/Ho ratios of 
30.5, 37.1, 33.9 and 36.4 for the banded shales, Fe-rich carbonates, Mn- 
rich sandstones, and Mn-rich black shales, respectively, and 36.3 for the 
two pyritized units (Fig. S2A & Table S3). The Fe-carbonates have the 
lowest EREE + Y content, moderate for the Mn-rich lithologies and 
highest for the banded shales (Fig. S2B & Table S3). These observations 
highlight a relational increase in total ZRFE + Y with the greater detrital 
composition of the banded shales and the higher Y/Ho ratios in the Fe- 
Mn-rich and pyritic sediments in association with a predominantly 
chemical origin (Fig. S2). The studied section overall shows a typical 
seawater-like positive shale normalized Y anomaly (Fig. 7D-G). 

However, the REE data, together with those of contemporaneous 
northwestern Congo Craton BIFs, associate the lithologies with mainly 
<250 °C low temperature hydrothermal fluids (Fig. 7C-G). Being a 
strong component of low temperature shallow-submarine hydrothermal 
fluids (Breuer and Pichler, 2013), arsenic is also enriched by several 
orders of magnitude in the Fe-rich carbonates compared to the under- 
lying banded shales with respect to PAAS (Fig. 3P), with normalization 
to detrital Al suggesting a authigenic or chemical precipitation (Fig. 41). 
Authigenic and or chemical enrichment of Eu and Ce is supported by a 
lack of meaningful correlation between detrital Al, Eu/Eu* and Ce/Ce* 
anomalies (Fig. S3A-B). With respect to PAAS composition, Eu and Ce 
are both several times enriched in the Fe-rich carbonates, with similar 
but muted enrichment in the Mn-rich black shales, sandstones and 
pyritized black shales relative to the detritus-rich banded shales 
(Fig. S3C-D). This observation further suggests a more obvious authi- 
genic and or chemical origin for the REEs in the upper metalliferous 
lithologies, relative to major detrital contribution in the underlying 
banded shales. 


4.3. Paleo-weathering depositional dynamics 


An average CIA value of 71.6 % for the studied section (Figs. 4K & 
5A), evocative of the 70-75 % UCC shale estimates, indicates overall 
moderate weathering of source rocks (Fedo et al., 1995; Fedo and 
Babechuk, 2023), being generally higher for the upper metalliferous 
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Fig. 3. Major and trace element enrichment patterns up the studied section. Enrichment factors (EF) were calculated according to Tribovillard et al. (2006). 
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Fig. 4. Authigenic elemental enrichment with respect to the Chemical index of alteration (CIA) weathering and detrital Al proxies. A-J) Trace metal(loid)s and 
phosphorus concentrations normalized to the Al detrital proxy plotted against depth (sequence/lithostratigraphy) and CIA. K) CIA trends across sequence/ 
lithostratigraphy. 
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Fig. 5. Bar and whisker diagrams showing the relationship between average CIA and Sequence/lithostratigraphy and total Fe (Fe203), Mn MnO, TiO and Al203 in 
the banded shales, Fe-rich carbonates and Mn-rich black shale units. A) Average distribution of measured CIA weathering indices along sequence/lithostratigraphy, 
from the underlying banded shales passing upwards through the thin pyritized layer, the Fe-rich carbonates and a second thin pyritized unit, capped by Mn-rich 
sandstones and black shales. Bars represent standard deviation from the mean. B) Relationship between TiO% content and the banded shales, Fe-rich carbonates 
and Mn-rich black shale units. C) Relationship between Al2O3 content and the banded shales, Fe-rich carbonates and Mn-rich black shale units. D) Relationship 
between Fe20; content and the banded shales, Fe-rich carbonates and Mn-rich black shale units. E) Relationship between MnO content and the banded shales, Fe-rich 
carbonates and Mn-rich black shale units. F) Relationship between Fe203/MnO ratios and the banded shales, Fe-rich carbonates and Mn-rich black shale units. G) 
Relationship between Fe203/TiO» ratios and the banded shales, Fe-rich carbonates and Mn-rich black shale units. H) Relationship between Fe203/Al203 ratios and 
the banded shales, Fe-rich carbonates and Mn-rich black shale units. I) Relationship between MnO/TiOz ratios and the banded shales, Fe-rich carbonates and Mn-rich 
black shale units. J) Relationship between MnO/AI,03 ratios and the banded shales, Fe-rich carbonates and Mn-rich black shale units. Bd-SH represents banded 
shales; Py-BS, Pyritized black shales; Fe-Carb, Fe-rich carbonates; Mn-SS, Mn-rich sandstones; Mn-BS, Mn-rich black shales. 


facies where CIA averaged 82.5 % for the Mn-rich black shales (Fig. 5A). term recycling and sorting. Trace element normalization to the Al 
As already mentioned, TE and ME enrichment correlate well with detrital proxy associates enhanced chemical weathering of source rocks 
sequence stratigraphic CIA records. The CIA data tend to plot along the with the deposition of the upper metalliferous facies (Fig. 4) and 
A-CN axis, parallel to expected weathering trends for UCC, TTG and increasing dissolution of various trace metals and P in seawater by 
granite, without significant deviation towards the A-K axis, suggesting several orders of magnitude (Fig. 3). 


limited post-depositional K enrichment (Fig. 8A). The samples cluster 
predominantly in the felsic field, mainly around UCC and TTGs on the 
compositional discriminant diagrams of Th/Cr versus La/C and Th/Sc 
versus La/Sc, with associated Al/Ti ratios signifying major contributions 
from felsic source rocks (Fig. 8B-E). A few samples occupy the limit 
between mafic and felsic crust, hinting on possible minor involvement of 
intermediate rock sources. Moreover, Th/Sc versus Zr/Sc and Al203 
versus SiO2 plots (Fig. S4), together with a mainly felsic crust prove- 
nance of the samples (Fig. 8), indicate that the CIA values and the 
geochemical profiles are not biased by sediment recycling, sorting and 
lithology above the banded shale boundary. On the other hand, the 
strong detrital composition of the banded shales is consistent with long- 


4.4. Redox depositional conditions 


A combination of Corg, Ccarbonates, N, true Ce anomalies, enrichment 
of redox sensitive TEs and Fe speciation systematics (Figs. 3, 9A-K, S5 & 
Table S4), delineate two broad redox trends across the studied section. 
Low and invariant bulk total Corg and N systematics evolved from the 
banded shale to being markedly enriched in the overlying Fe-Mn-rich 
facies (Fig. 9C, D). This boundary is associated with stable N isotopes 
switching towards more positive values and Corg isotopes towards more 
negative ratios (Fig. 9B, E). Inferred redox changes and Corg isotope 
distribution correlate well with those of the FBip., unit of the adjacent 
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Fig. 6. Cross plots showing the behaviour of Fe203/TiO2 and MnO/TiO, ratios 
with respect to Fe203 and MnO concentrations in the banded shales, Fe-rich 
carbonates and Mn-rich black shale units. A) Fe203/TiO2 ratios showing cor- 
relations for the individual banded shales, Fe-rich carbonates and Mn-rich black 
shale units. B) Fez03/TiO2 ratios showing correlations across the banded shales, 
Fe-rich carbonates and Mn-rich black shale units. C) MnO/TiO, ratios showing 
correlations for the individual banded shales, Fe-rich carbonates and Mn-rich 
black shale units. D) MnO/TiO2 ratios showing correlation across the banded 
shales, Fe-rich carbonates and Mn-rich black shale units. 


42 km? Bangombe Plateau (Canfield et al., 2013; Ossa Ossa et al., 2018). 
The banded shales associate with oxic bottom waters and the overlying 
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characterized by anoxic bottom waters experiencing ferruginous and 
euxinic conditions (Fig. 9G, H & S5). PAAS normalized REE patterns 
with negative Ce anomalies are muted in all lithologies, although 
varying levels of light REE (LREE) depletion and slightly positive La 
anomalies are discernible (Fig. 7D-G). Association of true Ce anomalies 
to oxic bottom oxic water deposition of the banded shales (Fig. 9F, K) 
coincide with oxygenated conditions inferred by Fe speciation analysis 
(Fig. 9G, H). Enrichment of a multitude of redox sensitive metal(loid)s 
relative to PAAS is inferred for the upper banded shale transition 
(Fig. 3H-J), with Mo, U and TOC enrichment suggesting hydrological 
restriction of the basin (Figs. 3 & 9I, J). 


5. Discussion 
5.1. Provenance and facies depositional mechanisms 


5.1.1. Weathering influence 

According to Fedo and Babechuk (2023), mafic igneous rocks rich in 
biotite, olivine, pyroxene and amphibole can influence the CIA index. 
However, provenance analysis suggests a mainly UCC-felsic-TTG affinity 
for our samples (Fig. 8). With average TTGs containing <5 wt% ferro- 
magnesian elements combined, implies such low ferromagnesian con- 
tent could not have skewed our observations. 

The flattening out of the positive correlation between Zr/Sc and Th/ 
Sc ratios in the banded shales (Fig. S4A) and the <35 Zr/Sc ratios in the 
upper metalliferous units, together with coeval high CIA values (Fig. 4K 
& 5A), indicate a primary composition controlled by similar source rock 
composition and negligible primary sediment reworking (McLennan 
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Fig. 7. Link between lithology and hydrothermalism inflence. A) Fe/Ti versus Al/Al + Fe + Mn hydrothermal proxy cross plot. Distinct hydrothermal influence is 
evident when Fe/Ti is >100, combined with low Al/(Al + Fe + Mn) ratios (Bostrom, 1970; Bostrom et al., 1972, 1973; Pecoits et al., 2009). B) Fe/Ti versus Al/Al + 
Fe + Mn plots for the Gouap BIFs of the northwestern margin of the Congo Craton (Soh Tamehe et al., 2021). Ti concentrations of the Mewongo BIFs were below 
detection limit to allow a similar comparison. C) Chrondrite normalized (CN) Sm/Y versus Eu/Sm cross plots (Alexander et al., 2008; Deassou Sezine et al., 2022). D- 
G) Rare Earth Element (REE) distribution compared with high and low temperature hydrothermal fluids, modern seawater, and the Paleoproterozoic BIFs of a similar 
age in the northwest margin of the Congo Craton (Soh Tamehe et al., 2021). Sm/Yb vs Eu/Sm ratios according to Alexander et al., 2008(53). Average hydrothermal 
fluid REE composition reconstructed according to Bau and Dulski (1999) and compared to modern seawater (Alexander et al., 2008 and references therein). Although 
the northwest Congo Craton BIFs have experienced metamorphism, their seawater REE and hydrothermal composition has been shown to be preserved (Soh Tamehe 
et al., 2021, 2022; Deassou Sezine et al., 2022). BS, black shale. SS, sandstone. CB, carbonates. SH, shale. 
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Fig. 8. Paleo-weathering (CIA) and Source provenance reconstruction of the studied facies. A) A-CN-K (Al203—CaO + Na2O-K20) ternary diagram (Fedo et al., 1995; 
Fedo and Babechuk, 2023). Arrows 1, 2 and 3 indicate potential weathering trends for basalt and andesite, UCC-TTG and granite, respectively. Data for basalt, 
andesite, TTG and granite are from Condie (1993). UCC and PAAS data are from Taylor and McLennan (1985). B) Th/Sc versus La/Sc provenance biplots. C) Th/Cr 
versus La/Cr provenance biplots. D) Th versus Sc provenance biplots. E) Al/Ti ratios plotted against sample depth and stratigraphy, showing that the entire set of 
samples have a strong terrigenous composition. According to this model, Al/Ti ratios closer to 20 reflect average continental rock values and those nearer 5, 
weathering products of oceanic rocks (Bostrom et al., 1969; Bostrom, 1970). 
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Fig. 9. Redox reconstruction and distribution of organic carbon, bulk nitrogen and their isotopic compositions in the studied cross section. A) Distribution of carbon 
isotopes in organic carbon and carbonates associated with Fe-Mn-rich carbonate deposits in the FB, Formation Member according to Weber and Gauthier-Lafaye 
(2013). B-C) Organic carbon systematics for the four drill holes analyzed. Inset in panel c shows the spread of the data excluding the pyritized sample with total 
organic carbon content of up to 31 wt%. D-E) Nitrogen systematics for the four drill holes studied. F) Ce/Ce* versus Pr/Pr* anomalies calculated according to Bau and 
Dulski (1996). G-H) Iron speciation across sequence stratigraphy. I) Mo versus U enrichment patterns for the four drill core holes studied (Tribovillard et al., 2006; 
Algeo and Tribovillard, 2009). J) Mo vs TOC cross plot for the drill core holes studied compared to modern marine sediment profiles with different redox and basin 
restriction properties (Algeo and Lyons, 2006). K) Ce/Ce* versus Pr/Pr* plot for the four drill core holes studied. 


accompanied by enrichment of poorly mobile heavy minerals (notably clay-size minerals rich in Al2O3. The end point of this process can be 
Zr in zircon minerals), while accumulation of stable minerals like quartz expressed as rising Si02/A1,03 ratios coinciding with falling TE content 
dilutes sediment components and elements (McLennan et al., 1993). in facies affected by extensive sediment reworking. Although we see 
Immobile Th and Sc on the other hand do not undergo significant systematic enrichment of SiOz and Al203 progressing from the banded 
fractionation during transport and deposition (McLennan et al., 1993). shales to the upper metalliferous sediments (Fig. S4B), the Si02/Al,03 

According to McLennan et al. (1993), the >35 Zr/Sc ratios in the ratios suggest SiOz is more enriched in the Fe-rich carbonates than Al203 
banded shales are more consistent with primary sediment sorting and (Fig. S3C). However, the elevated TE content of the Fe-rich carbonates 
recycling, i.e., reworking which enriches quartz (SiO2) at the expense of (Figs. 3-4) does not support primary sediment reworking as the source of 
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SiO» enrichment in this section. Low Si02/Al203 ratios in the Mn-rich 
lithologies correspond well with elevated TE composition and negli- 
gible sediment reworking as corroborated by the Th/Sc versus Zr/Sc 
cross plots (Fig. S4A). Weak SiO2/Al203 ratios in the banded shales 
(Fig. S4C) due to co-enrichment of SiOz and Al2O3 (Fig. S4B) are 
consistent with low TE content (Fig. 3) attributed to primary sediment 
reworking (Fig. S4A). Overall, succession of low CIA values in the 
banded shales by stepwise increase up section, suggest the CIA values 
and geochemical signals of the metalliferous layers reflect increasing 
chemical weathering more than change in source rocks and primary 
sediment reworking. 


5.1.2. Hydrothermal influence 

Deposition of the Fe-Mn-rich rocks has been previously attributed to 
upwelling of deep-sea hydrothermal Fe?+/Mn**-rich fluids to conti- 
nental margin waters (Gauthier-Lafaye and Weber, 2003; Ossa Ossa 
et al., 2018). However, contribution by continental weathering and the 
source of hydrothermal fluids remain to be co-elucidated. Therefore, we 
use the Fe/Ti versus Al/(Al + Fe + Mn) ratio to distinguish the primary 
influence of hydrothermal activity from other sources. This assumption 
is reliant on the premise that calibration of this proxy against average 
continental and oceanic crust, active hydrothermal ridge sediments and 
a mixture of these sources, associates Fe/Ti > 100 and low Al/(Al + Fe 
+ Mn) ratios with hydrothermal fluid influence (Bostrom, 1970; 
Bostrom et al., 1972, 1973; Pecoits et al., 2009). Accordingly, Al/(Al + 
Fe + Mn) ratios correlate geography with submarine volcanic ridge 
activity and the distribution of pelagic seawater Al/(Al + Fe + Mn) ra- 
tios (Bostrom et al., 1969; Bostrom, 1970; Bender et al., 1971). 

Our >100 Fe/Ti ratios appear to suggest stronger involvement of 
hydrothermal activity in the deposition of the Fe-rich carbonates than in 
the later stage formation of the Mn black shales, with little or no 
contribution indicated for the banded shales. However, Mn black shale 
and banded shale Fe/Ti ratios are more diluted by Ti loading compared 
to the Fe-rich carbonates (Fig. 6A), because of greater contamination by 
continental detritus. Nevertheless, the Mn black shale Al/(Al + Fe + Mn) 
ratios are distinctively lower than for the detritus-richer banded shales. 
The most parsimonious interpretation of these trends is that both the 
Mn-rich black shales and Fe-rich carbonates experienced more influen- 
tial hydrothermal activity than the banded shales. However, deposition 
of the Mn-rich black shales appears to have coincided with gradual 
waning of hydrothermal activity and increasing detrital supply. This 
association is supported by initiation of Mn enrichment in sandstones 
underlining the Mn black shales. Geochemically, the Fe-rich carbonates 
contain maximum Al/(Al + Fe + Mn) ratios of ~0.2, doubling to >0.5 in 
the Mn black shales and up to >0.8 in the banded shales. Moreover, 
connection of the Fe-rich carbonates with greater hydrothermal fluid 
influence is supported by (Sm/Yb)cn versus (Eu/Sm)cn cross plots in 
which up to ~1-5 % hydrothermal fluids may have been involved, with 
a maximum of ~0.1-1 % for the Mn black shales (Fig. 7C). 

To further constrain the source and potential hydrothermal fluid 
temperatures, we consider the estimated depositional age of 
~2100-2000 Ma for the FB,, unit, which is believed to coincide with 
collision of the Congo-Sao Francisco cratons during the Eburnean- 
Transamazonian orogeny (Pambo et al., 2006; Barbosa and Sabaté, 
2004; Teixeira et al., 2017; dos Santos et al., 2023). Pambo et al (2006) 
describe a model which juxtaposes the Congo Craton to the Sao Fran- 
cisco and the Transamazonian Cratons in Brazil, related to a conver- 
gence activity and volcanism, contemporaneous with the onset of 
orogenesis at ~2100-2040 million years ago. In this model, they pro- 
pose a broad lack of relief in the Francevillian basin due its location over 
the subducted crust and therefore out of the orogenic range. Compari- 
sons with modern high and low temperature hydrothermal fluid REE 
regimes indicate that this tectonic event may have resulted in the 
deposition of low temperature hydrothermal BIFs on the northwestern 
Congo Craton margin, Northwest of the Gabonian block to which the 
Francevillian sub-basin belongs (Soh Tamehe et al., 2021, 2022; 
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Deassou Sezine et al., 2022). The ~2100 Ma Gouap BIFs are thought to 
represent potential deeper water equivalents of the shallow continental 
margin Francevillian Mn-rich facies (Soh Deassou Sezine et al., 2022; 
Soh Tamehe et al., 2022). Our provenance analysis however affiliates 
the studied section to UCC and TTGs rocks (Fig. 8), although they also 
display low temperature hydrothermal axial-like fluid Ca enrichment 
signal and overall relative depletion of K and Mg (Fig. S6). 

Compositionally, the northwestern Congo Craton BIFs and the Fe- 
rich carbonates are moderately depleted in LREEs relative to heavy 
REEs (HREEs) and express positive Eu anomalies (Fig. 7E). The precip- 
itation of Eu from >250 °C hydrothermal fluids (Bau and Dulski, 1999; 
Alexander et al., 2008; Soh Tamehe et al., 2021, 2022; Deassou Sezine 
et al., 2022), imprints positive Eu anomalies when high temperature 
hydrothermal fluid contribution is >0.035 % (Delvigne et al., 2012). 
However, our broadly mild Eu anomalies are more akin to the gener- 
alized records of Proterozoic seawater trends (Planavsky et al., 2010). 
Further, the true Ce anomalies reflect those of typical Archean- 
Paleoproterozoic Fe formations (Planavsky et al., 2010). Our true posi- 
tive Ce anomalies and Fe speciation data are consistent with previous 
suggestions that the reworked underlying banded shales formed beneath 
full to partially oxygenated bottom waters, and the succeeding Fe-rich 
carbonates and the Mn-rich black shales and sandstones beneath pre- 
dominantly anoxic (ferruginous/euxinic) waters (Canfield et al., 2013; 
Ossa Ossa et al., 2018). 

On the other hand, the northwestern Congo Craton Mewongo BIFs 
display negative Y anomalies absent in the Gouap BIFs, while positive 
seawater-like Y anomalies are discernible in the Francevillian facies. The 
BIFs fall mainly between (Sm/Yb)cn and (Eu/Sm)cn fields indicating 
dominant seawater influence and ca. <0.1 % hydrothermal fluid 
contribution (Soh Tamehe et al., 2021, 2022; Deassou Sezine et al., 
2022). The Francevillian Fe-Mn-rich facies show a similar dominant 
seawater signal, with the Fe-rich carbonates being more inclined to 
greater hydrothermal fluid influence than the Mewongo and Gouap BIFs 
(Fig. 7C). However, REE compositional analysis suggests that the Gouap 
BIFs and the Francevillian lithologies may be more identical than the 
Mewongo BIFs (Fig. 7D-G). Out of the three Paleoproterozoic BIFs 
described in the northwestern Congo Craton, all containing significantly 
lower P205 concentrations compared to the Francevillian Fe-rich car- 
bonates, the Gouap BIFs exhibit the highest bulk P205 content (Fig. S7). 
Association of these deposits with low temperature hydrothermal fluids 
probably reflects seawater mixing with arc-related hydrothermal fluids, 
following westward subduction of the Sao Francisco Craton underneath 
the West Gabonian block during the Eburnean-Transamazonian orogeny 
(Soh Tamehe et al., 2022). 


5.2. Basin architecture influence 


The CIA values and inferred correlation to upper UCC and TTG 
source rock composition, combined with exceptional P enrichment in 
the Fe-rich carbonates (Figs. 3M & 4I), agree with the exposure of 
3050-2700 Ma felsic TTG-like granitoids in the northern and southern 
parts of Gabon, considered to be protoliths of the Francevillian Group 
facies (Bankole et al., 2020). These rocks share similar Archean ages of 
3200-2900 Ma with TTG suites in northwestern Congo and the Sao 
Francisco Craton (Barbosa and Sabaté, 2004; Alexander et al., 2008; 
Deassou Sezine et al., 2022; Teixeira et al., 2017). 

Mo, U and TOC enrichment patterns indicate a hydrologically 
restricted basin architecture for the Francevillian sub-basin (Fig. 9I, J), 
in agreement with suggestions by Reynaud et al. (2018) and Aubineau 
et al. (2021), which contrasts with the foreland open ocean facing Nyong 
Complex of the northwestern Congo Craton. For example, the increasing 
enrichment of Corg at the expense of Mo (Fig. 9J), especially in the 
anoxic Mn-rich black shales and Mn-rich sandstones, is indicative of 
strong removal of Mo from a redox stratified water column, without 
replenishment from deep water sources (Algeo and Lyons, 2006). Mo- 
lybdenum and Corg are well correlated in the underlying oxic banded 
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shales (R? = 0.7490), followed by no correlation in the overlying Mn- 
rich Fe carbonates (R? = 0.0052), before transitioning into a negative 
linear correlation in the Mn-rich black shales (R? = 0.7002) (Fig. S8). 
These trends are consistent with a change from oxic bottom water 
conditions towards more redox stratified waters reminiscent of an 
epicontinental sea setting like the Black Sea (Algeo and Lyons, 2006; 
Algeo and Tribovillard, 2009; Dong et al., 2023a,b). The positive cor- 
relation between Mo and Corg in the pyritized black shales possibly re- 
flects pyrite’s high affinity for Mo in the more intense anoxic conditions 
that enhanced co-preservation of Corg, rather than an increase in 
seawater Mo concentrations. 

Furthermore, the Francevillian sub-basin and the open ocean facing 
Nyong Complex continental platforms show distinct redox attributes, 
unique hydrothermal dynamics and different degrees of metamorphic 
and/or diagenetic influences, suggestive of distinct paleogeographic 
settings that can explain their divergent elemental enrichment patterns. 
For example, P205 content in the Gouap and Mewongo BIFs average 
0.18 and 0.06 wt%, respectively (Soh Tamehe et al., 2021; Deassou 
Sezine et al., 2022), 1.35 wt% for the Francevillian Fe-rich carbonates 
and 0.3 wt% for the Mn-rich black shales (Table $3). Further, the Fe-rich 
carbonates contain a lower average Fe content of 33.6 wt% compared to 
>41 wt% for the Gouap and Mewongo BIFs (Soh Tamehe et al., 2021; 
Deassou Sezine et al., 2022). Similarly, the Gouap and Mewongo BIFs 
record average MnO content of 0.15 wt% and 0.03 wt%, respectively, 
relative to up to 6.1 wt% for the Fe-rich carbonates and 23.7 wt% for the 
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Mn black shales and sandstones. Moreover, the northwestern Congo 
margin craton BIFs were deposited in relatively oxic deeper waters (Soh 
Tamehe et al., 2021; Deassou Sezine et al., 2022), while the Fe-Mn-rich 
Francevillian facies formed in shallow continental margin anoxic 
waters. 

These paleogeographic and tectonic setting influences imply that 
Francevillian sub-basin and the northwestern Congo Craton were likely 
connected to different submarine hydrothermal fluid sources and local 
geological and weathering processes, both controlled by the onset of the 
Eburnean-Transamazonian orogen. 


5.3. Diagenetic influence 


Multiple lines of evidence, including negative 5'°Ccarbonate Values in 
the Mn-rich black shales and siderite in the Fe-rich carbonates, a switch 
to Corg enrichment (Fig. 9A, B) and coupling to more negative ois 
(Fig. 9B) and more positive 5'°N values (Fig. 9E) and Ce anomalies 
(Fig. 9F), anoxia and euxinia (Fig. 9G, H), collectively point to early 
diagenetic conditions involved in biotic and abiotic alteration of a pre- 
cursor Fe-Mn oxide reservoir to Fe-Mn carbonates. For example, low 
negative 5'?Ccarbonate Values of —10 %o for the Mn carbonates and —13.4 
%o for siderite in the Fe-rich carbonates (Fig. 9A), largely suggest a 
coupling of early diagenetic microbial reduction of Fe-Mn oxides, sulfate 
and nitrate to the oxidation of Corg (Fig. 10B). This is in addition to 
previously postulated methanotrophic methane oxidation (Weber and 
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Fig. 10. A conceptual model showing the link between Fe and Mn cycling, redox conditions and hydrothermal activity in the Francevillian sub-basin. A) Fe and Mn 
cycling in a fully oxic shallow water basin before hydrothermal activity leading to long-term scavenging of P, trace elements (TEs) and rare earth elements (REEs) 
from seawater by Mn oxides and Fe oxides. B) Water column redox stratification induced by strong hydrothermal activity leading to reductive weathering and 
recycling of the P-TEs-REEs-rich Fe and Mn oxides, back to seawater. Further dissolution of the Fe-Mn deposits, facilitated by sulfide produced by the sulfate reducing 
prokaryotes, increased dissolved seawater sulfide concentrations, releasing more TEs and P to seawater. This would have been augmented by continental weathering 
and TEs supply by hydrothermal emissions. The continued burial of Fe and Mn as carbonates would have reduced the ability of the Fe and Mn oxides to recapture 


seawater P and bioessential TEs, thereby increasing bioavailability. 
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Gauthier-Lafaye, 2013). The decomposition of microbial mats in modern 
settings depletes 13¢ in carbonates (Sanz-Montero et al., 2023). More- 
over, the Mn carbonate and siderite O Carbonate fall within the range of 
diagenetic biogenic methane production and Corg oxidation in modern 
ferruginous environments (Vuillemin et al., 2023). 

Reduction of Fe oxides by the dissimilatory Fe-reducing bacteria 
with Corg under anoxic conditions produces biogenic 8!3C-depleted 
siderite (FeCO3) with isotopic values of <—10 %o (Zhang et al., 2001). 
Low FeCO3 6!°C values could also originate from dissolved inorganic C 
(DIC) in hydrothermal-seawater-mantle fluid mixtures with a 5'3C of ca. 
—6.5 %o or from kinetic isotope fractionation when FeCOs3 precipitates 
directly from seawater-derived DIC (Jiang et al., 2022). For instance, the 
observed shift to positive Ce/Ce* anomalies in the Fe-Mn-rich lithol- 
ogies, is a diagenetic MnCO3 and FeCO3 enrichment signature (Huang 
et al., 2022; Jiang et al., 2022; Dong et al., 2023a,b; Mhlanga et al., 
2023). Typically, Ce prevails as tetravalent ions in oxygenated seawater, 
being less soluble and more reactive with Mn and Fe oxides than other 
REEs (Bau and Dulski, 1999). The ensuing <1 Ce/Ce* in oxygenated 
seawater is reflected in the oxic banded shales, being >1 in the anoxic 
metalliferous facies. In this scenario, highly reactive Fe (Feyp) to total Fe 
(Fer) and pyrite Fe (Fepy) to Feyr ratios mirror comparable bottom water 
anoxia for the correlative Bangombe Plateau (Canfield et al., 2013). 

A similar diagenetic model involving the reduction of a seafloor 
ferromanganese oxides linked to changing redox conditions, was 
recently proposed for dolostones associated with Fe-Mn enrichment in 
the correlative FB Formation in the Lastoursville sub-basin (Yoshida 
et al., 2023). This observation suggests a basin-wide mechanism for the 
deposition of Fe and Mn in the Francevillian group. Differences in 
mineralogical enrichment patterns between the Lastoursvillian and 
Francevillian sub-basins are probably related to a complicated basin 
architecture that resulted in the deposition of the Francevillian Group in 
four distinct rift basins probably characterized by horsts and grabens, in 
Franceville, Okondja, Lastoursville and Booué (Weber, 1968; Feybesse 
et al., 1998; Pambo et al., 2006; Yoshida et al., 2023). 

In this diagenetic model, formation of the precursor Fe-Mn oxides 
(Fig. 10A) was probably aligned with the widespread oxic conditions 
that preceded deposition of the Fe-rich carbonates (Yoshida et al., 
2023). Moreover, the observed primary sediment reworking that marks 
the banded shales implies progressive loss of soluble elements into so- 
lution could have contributed to systematic enrichment in seawater and 
subsequent incorporation and precipitation into precursor seafloor Fe- 
Mn oxides. The Fe-Mn oxides were then later transformed to Fe-Mn 
carbonates when anoxic hydrothermal fluids overwhelmed bottom wa- 
ters and sea level rise (Canfield et al., 2013; Ossa Ossa et al., 2018; 
Reynaud et al., 2018; Yoshida et al., 2023) promoted shoaling of anoxic 
conditions (Fig. 10B). In addition to this already existing seafloor Fe-Mn 
deposit, the onset of hydrothermal activity likely supplied additional 
hydrothermal Fe and Mn that was oxidized in the redoxcline to form 
surplus Fe-Mn oxides that was further transformed into Fe-Mn carbon- 
ates beneath the redoxcline and in the unconsolidated sediment pile 
during early diagenesis (Fig. 10B). It is well-established that Fe-Mn 
oxide reductive processes in unconsolidated sediments enrich P, Mn?* 
and Fe** in overlying anoxic bottom waters. 

The seafloor altered Fe-Mn sources were likely separated both 
temporally and spatially, considering that under anoxic conditions, 
faster reductive dissolution of Mn over Fe oxides should first result in 
preferential deposition of the Mn-rich black shales and sandstones 
before the Fe-rich carbonates, which is not observed. It is also possible 
that formation of some of the Fe-Mn carbonates could have bypassed the 
Fe-Mn oxide phase if alkalinity produced by microbial oxidation of Corg 
beneath the redoxcline was sufficient to allow this to happen (Roy, 
1992; Wittkop et al., 2020). 

Our data suggest that chemical weathering of an inferred felsic 
protolith proximal to or in shallow submerged oxygenated platforms, 
probably allowed oxidative weathering to progressively enrich the 
seafloor with Fe-Mn oxides and co-drawdown of P and bioessential TEs 
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and REEs from seawater (Fig. 10A; Klinkhammer et al., 1983; Calvert 
and Pedersen, 1996; Wheat et al., 1996, 2003; Feely et al., 1998; Kuhn 
et al., 1998; Yoshida et al., 2023). Under these oxic conditions, V, U and 
Mo would have been released by oxidative weathering into seawater 
(Tribovillard et al., 2006). In conditions accompanied by strong oxida- 
tion of reduced Fe and Mn and high Corg content, Mo, V and U are co- 
enriched in Fe-Mn oxide/Co:,-rich sediments (e.g., Anderson, 1982; 
Hein et al., 1997). Increasing enrichment of Mo over U is expected in 
bottom water euxinic sediments (Fig. S5; Algeo and Maynard, 2004; 
Tribovillard et al., 2006). 

Further, the high affinity of Ce(IV) formed in oxic waters to bind Fe- 
Mn oxides, Corg and clays, would have resulted in shuttling to bottom 
waters and eventual burial in seafloor sediments (Byrne and Sholkovitz, 
1996; Kuhn et al., 1998). This interpretation is supported by co- 
enrichment of Corg, REEs and various redox sensitive trace metals in 
the metalliferous horizon. It is thought that remobilization from similar 
sources concentrated REEs in Paleoproterozoic redox stratified bottom 
waters (Planavsky et al., 2010). Like U, V and Mo enrichment, positive 
Ce anomalies characterize the upper banded shale boundary that marks 
a sharp transitional boundary from oxic to redox stratified water mass. 
Moreover, the recognized absence of positive Ce anomalies in hydro- 
thermal fluids (Klinkhammer et al., 1983; Kuhn et al., 1998) supports 
the shuttling of Ce by potential Fe-Mn oxides to the seafloor, followed by 
later diagenetic enrichment in the Fe-Mn carbonates. 

The possibility also exists that the Mn-rich black shales formed 
through partial ventilation and alteration of an unidentified sulfidic 
sediment source according to the model of Dong et al. (2023a,b). 
However, large scale Mn carbonate deposits are predominantly formed 
in anoxic sediments below oxygenated surface waters in nearshore 
continental margin environments, where they are enriched via a surface 
Mn oxyhydroxide shuttle (Calvert and Pedersen, 1996). Moreover, 
sequential deposition of the Fe-rich carbonates in combination with 
greater seawater nutrient enrichment is more agreeable with immediate 
deposition of Corg-rich pyritized black shales above the Fe-rich carbon- 
ates. As explanation, the ensuing enrichment of seawater with P and 
trace metal nutrients would have stimulated intense primary production 
in surface waters (Fig. 10B). Reduced hydrothermal fluid intrusion, 
together with high surface water Corg production would have intensified 
further deoxygenation of deeper waters by coupling microbial con- 
sumption of oxygen to the oxidation of Corg, in turn enabling microbial 
pyrite deposition. According to this model, precipitation of the Mn-rich 
facies would only have become permissible when these extreme 
reducing conditions abated to enable deposition. 


5.4. Nutrient enrichment influence 


Phosphorus could have been supplied by hydrothermal plumes 
following alteration of seafloor basalts by anoxic hydrothermal fluids 
(Syverson et al., 2021; Rasmussen et al., 2021), in addition to major 
riverine delivery to surface waters (Sharoni and Halevy, 2022). How- 
ever, the probability of enriching up to ~4 wt% P in the Fe-rich car- 
bonates by alteration of seafloor basaltic rocks by anoxic hydrothermal 
fluids alone—i.e., bypassing the Fe oxide enrichment phase—is not 
supported by our provenance analysis. Moreover, it is suggested that 
seafloor weathering on its own contributes only a small fraction of P to 
seawater with the chemical weathering of continental silicates being the 
dominant input (Wheat et al., 1996, 2003, 2008; Sharoni and Halevy, 
2022). The latter observation aligns with lack of evidence pointing to 
major P supply to paleo-marine waters by seafloor weathering to a level 
comparable with the Francevillian records (Fig. 11; Reinhard et al., 
2017). Otherwise, this mechanism which would have been more prev- 
alent in anoxic Archean-Paleoproterozoic seawater subjected to greater 
hydrothermal fluid influence, should have recorded similarly high P 
levels in marine sedimentary records spanning the Precambrian oceans 
(Planavsky et al., 2010). 

Moreover, the greater chemical weathering indices recorded in the 
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Fig. 11. Log normal comparison of Zn and P composition for samples with published marine shale values through Earth history. A) Zn concentration in ppm. B) Zn 
concentration normalized to Al as a detrital and authigenic enrichment proxy. C) Sample P data plotted in weight percent P concentrations compared to marine shales 
through Earth history. The Zn and P data are adapted from Scott et al. (2013) and Reinhard et al. (2017), respectively. 


Mn black shales and Fe-rich carbonates are not coeval with comparable 
P205 enrichment despite their similar rock source provenance. There- 
fore, a change in source rock composition cannot explain the distinct 
coupled Fe, Mn and P enrichment patterns. Likewise, a mainly detrital 
origin for apatite concentrated in the Fe-rich carbonates (Mossmann 
et al., 2005) cannot be justified since the Mn black shales and sandstones 
show a greater detrital influence than the Fe-rich carbonates. Moreover, 
apatite enrichment in the Mn-rich black shales should be greater if the 
more negative 5 3Corg values in the Mn black shales and sandstones are 
indicative of greater microbial remineralization of P-rich organic mat- 
ter. However, this is not the case. When viewed in this light, enrichment 
of up to 4 wt% P in the Fe carbonates suggest unique feedback mecha- 
nisms existing between the P and Fe biogeochemistry, decoupled from 
Mn cycling. The absence of correlation between P and TOC concentra- 
tion (Fig. S9A) is consistent with low C/P ratios in the Fe-rich carbonates 
(Fig. S8B), likely reflects lower sedimentation rates and oxidation of Corg 
(Ingall and Van Cappellen, 1990). 

The most parsimonious mechanism for P enrichment in the Fe-rich 
carbonates therefore points to diagenetic alteration of a distinct sea- 
floor P-rich Fe oxide reservoir by dissimilatory iron reducing microbes 
and sulfide produced by sulfate reducers in bottom water reducing 
conditions. In this model, systematic burial of Fe carbonates and pyrite 
retained P in solution at the expense of Fe oxide dissolution (Fig. 10B). 
The resulting saturation of dissolved P at the seawater-sediment inter- 
face would have triggered apatite precipitation (Ruttenberg and Berner, 
1993; Ruttenberg, 2003). 

Our bulk Zn, Cu and Mo concentrations are consistent with historical 
paleo-dynamic marine black shale archives (Figs. 11A & S10-S11). 
However, much lower Zn/Al ratios across the sampled section suggest 
that seawater dissolved Zn concentrations may have been much inferior 
in the Francevillian sub-basin compared to similar contemporaneous 
settings and that the upper banded shale boundary coincided with a rise 
in dissolved seawater Zn content (Fig. 11B). Despite this overall lower 
dissolved seawater Zn content, quantitative Zn enrichment and frac- 
tionation of light Zn isotopes is recorded in macrofossils located in the 
strata immediately above the studied facies (El Albani et al., 2023; Ossa 
Ossa et al., 2023). This observation appears to fall in line with the hy- 
pothesis that the evolution of biological Zn usage was perhaps controlled 
by cellular requirement than dynamic changes in seawater levels 
(Robbins et al., 2013). For instance, Zn biolimitation in some 
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cyanobacteria is managed by replacing Zn with Co (Saito et al., 2003), 
which tends to be more enriched relative to Zn at the upper banded shale 
boundary (Fig. S12). Similarly, some unicellular eukaryotes adapt to Zn 
deficiency by upregulating the expression of high-affinity Zn uptake 
systems and substituting Zn for Co and Cd (Timmermans et al., 2001; 
Kellogg et al., 2020, 2022). 

Bulk P enrichment in our samples is comparable to coeval records in 
this interval (Fig. 11C). The data highlight potential seawater P con- 
centrations equivalent to or nearing Ediacaran marine biosphere levels 
and beyond, a time when undisputed metazoans first appeared in the 
marine sedimentary rock record (Reinhard et al., 2017; Laakso et al., 
2020; Dodd et al., 2023). Considering the biogenic bridge between 
environmental change, seawater P content, trophic level complexity and 
body size with oxygenation (Laakso et al., 2020; Wang et al., 2020; Dodd 
et al., 2023), it is quite likely that localized environmental change in 
Franceville, together with a rise in dissolved seawater P content and 
oxygenation, might explain the subsequent explosion of aerobic 
macroscopic life in the overlying FB2 formation in the Francevillian sub- 
basin. 

The associated P/Al ratios point to maximum authigenic P enrich- 
ment during deposition of the Fe-rich carbonates, intermediate in the 
Mn-rich carbonates, sandstones and pyritized shales, and lowest in the 
banded shale interval (Fig. S13). The latter probably reflects contribu- 
tion by sediment reworking (Fig. S4A). Our fine-grained shales record 
sediment geochemistry in a near continental margin marine environ- 
ment, associated with significant authigenic P enrichment together with 
various bioessential metals in shallow marine waters. This contrasts 
with the relatively insignificant authigenic P enrichment reported by 
Reinhard et al. (2017) in shales of a similar age sampled from elsewhere 
(Fig. 11C), pointing to a unique local nutrient enrichment event in the 
Paleoproterozoic Francevillian sub-basin. The lower P content of the 
Mn-rich facies (Fig. 3M) further indicates that these seawater eutrophic 
conditions began to equilibrate up sequence to enable bottom water 
ventilation. The appearance of full oxic bottom waters in the immediate 
overlying FB2 Member (Canfield et al., 2013; Ossa Ossa et al., 2018), 
together with macrofossils (El Albani et al., 2010, 2014, 2019), and the 
proliferation of benthic microbial mats (Aubineau et al., 2018, 2019, 
2021), represent a potential hallmark of this transitionary boundary. 


E. Chi Fru et al. 
5.5. Nitrogen cycling influence 


Nitrogen and carbon systematic (Figs. S14-S15), showing a moderate 
linear covariation (R? = 0.35) between increasingly negative Caig and 
positive 5'°N enrichment (Fig. 9B, E & S14B), appear to suggest a weak 
but enhanced coupling between heterotrophic assimilation and reduc- 
tion of 4N-12C-rich Corg. This is because non-quantitative reduction of 
fixed !4N and elimination as N3 from bottom anoxic waters and sedi- 
ments, enrich sediments with !°N as a by-product. 

As a whole, oceanic planktonic biomass has an average molar C/N 
ratio of ~7 (Redfield, 1958), with biomass remineralization in the water 
column and sediments producing higher sedimentary molar C/N ratios. 
The predominantly 2 to 91 M C/N ratios in our samples (Table S6; 
Fig. S15) are consistent with those of other Francevillian facies (Aubi- 
neau et al., 2021; Kipp et al., 2018; Ossa Ossa et al., 2022), and hint at 
diagenetic reworking of primary biomass. In addition, remineralization 
of N-rich organic matter may release ammonium (NH) to pore waters 
that can be subsequently incorporated into phyllosilicates via substitu- 
tion for K cations, with further fractionation of d!°N (Williams et al., 
1995). Nonetheless, the lack of apparent correlation between TN and K 
(Fig. S14A), suggests negligible influence for phyllosilicate-bound N 
phases on the ð!?Npu values. 

The C and N isotope composition in sedimentary rocks is also sen- 
sitive to thermal alteration induced by diagenesis and metamorphism 
(Ader et al., 2006). However, the Francevillian sub-basin is well-known 
for little or no metamorphic overprinting. During thermal diagenesis, 
breakdown and volatilization of organic molecules enriches the residual 
pool with heavy C and N isotopes (Bebout and Fogel, 1992; Stiieken 
et al., 2017). This enrichment causes a shift in 0!°Npu values by <1 %o 
in lower greenschist facies that are not present in our samples. More- 
over, increasing 0!°Npurk value should coincide with decreasing 9 Corg, 
and it is expected that thermal alteration should promote N depletion 
relative to C as well as coeval decrease in TOC and !C. These conditions 
are typically represented by strong covariation between 0 Corg and 
TOC, 0 Corg and C/N, 2'°Nourk and TN, o’°’Nobutk and C/N, as well as 
0 Corg and d!°Nputk, none of which is observed across sampled section 
(Fig. S14B-F). The identification of randomly ordered diagenetic illite/ 
smectite mixed-layer clays (i.e., persistence of smectite) in <2 um clay- 
size fraction in the Mn black shales, have been interpreted as evidence 
for lack of significant post-depositional thermal alteration (Ngombi 
Pemba et al., 2014; Ikouanga et al., 2023). 

Modern oceans are dominated by aerobic N cycling in which nitrate 
is the dominant dissolved species (Stiieken et al., 2016). Nitrification 
followed by non-quantitative denitrification and/or anammox (anaer- 
obic ammonium oxidation coupled to nitrite reduction) in oxygen 
minimum zones are microbial N pathways invoked to explain modern 
positive '°N values (Dalsgaard et al., 2005; Sigman et al., 2009). On the 
other hand, biological Nə fixation preserves ð!®Npu values of —2.0 and 
+1.0 %o (Bauersachs et al., 2009; Zhang et al., 2014). High seawater 
ammonium (NH4) concentration in anoxic waters, together with intense 
uptake of NHj into biomass, can tune down biomass 5'°N towards more 
negative values (Wang et al., 2023). Nonetheless, it is well established 
that increasing positive 5'°N values correlate broadly with a large ni- 
trate pool coupled to active but non-quantitative denitrification and 
anoxia (Wang et al., 2023). The switch to more positive 5'°N coinciding 
with greater N enrichment in the metalliferous facies therefore suggests 
presence of a sizeable seawater nitrate reservoir due to enhancement of 
N fixation in surface waters. In this context, a rise in surface water pri- 
mary productivity, oxygenation and bottom water anoxia should be 
consistent with the shift towards sedimentary enrichment of N and more 
positive 5'°N (Fig. 9D, E). This observation agrees well with the two-step 
dynamic redox change documented by Fe, Ce and redox sensitive trace 
metal trends. 

Considering invariable atmospheric N isotope composition over the 
last ~3.0 billion years (Marty et al., 2013; Sano and Pillinger, 1990), 
positive 0'°Npuk Values in the range of +2.0 and +5.0 %o (average +3.2 
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+ 1.2 %o, n = 10) in the Francevillian FB, unit banded shales could 
reflect the preservation of an aerobic N cycle coupled to reductive N 
cycling. The Feyp/Fer ratios <0.22, together with low enrichment levels 
of redox-sensitive Mo, V, and U in the FB,, banded shales, support 
deposition under oxygenated bottom waters. Nonetheless, the 0 ’Nbulk 
values of the FB1» unit are lower compared to those varying between +4 
and +8 %o in the Francevillian FC Formation, linked to pervasive and 
stable aerobic N cycling (Kipp et al., 2018; Ossa Ossa et al., 2022). 

Bioavailable Fe, P, and N are primary biolimiting nutrients that 
control biological productivity over time (Capone and Hutchins, 2013; 
Reinhard et al., 2017; Tyrrell, 1999). We observe a moderate P enrich- 
ment (mean EFp = 1.7, n = 20) in the banded shales, while Fe is 
moderately depleted (mean EFpe = 0.5, n = 23) with respect to PAAS. 
This indicates that Fe limitation in the FB1» offshore environments could 
have hampered rates of marine productivity to some degree. Typically, 
N fixation by the nitrogenase enzyme is dependent on V, Fe and Mo 
(Hernandez et al., 2008; Zhang et al., 2014; Darnajoux et al., 2019). Our 
data show that these metals were enriched in seawater upon transition 
from the banded shales into the Fe-rich carbonates (Fig. 3B, H, J). 
Combined with the fundamental role of Mn in the splitting of water 
during oxygenic photosynthesis (Fischer et al., 2015) and the critical 
role of P in primary production (Tyrrell, 1999), the nutrient-rich con- 
ditions at the upper banded shale boundary would have intensified 
primary production in surface waters, intense N fixation, abundant Corg 
production and surface water oxygenation. These eutrophic conditions 
would have promoted the enrichment of positive but low 0° Npuk Values 
(<5%o) in the FB1» unit. This is because preferential nitrogen fixation by 
the Mo dependent nitrogenase enzymes does not strongly fractionate N, 
while strong denitrification of nitrate and anammox imparts a positive 
a°°N to residual sediment N content (Stiieken, 2013; Wang et al., 2023). 

The overlying FB; lithologies, marked by higher dNbutk values 
(average +5.7 + 0.9 %o, n = 21), together with ferruginous and sulfidic 
anoxic bottom water conditions resulting from upwelling of trace metal- 
P-rich fluids (Figs. 3-5), reflect an intensification of deoxygenation of the 
water column through partial loss of nitrate during denitrification at the 
redoxcline. The tremendous discharge of nutrients and anoxic waters 
would have likely shifted the redoxcline upwards. This basin-scale 
deoxygenation in the FB, unit would have subsequently led to the 
observed short-lived nitrate scarcity in the overlying FB Member, as 
hypothesized by Aubineau et al. (2019) and confirmed by Ossa Ossa 
et al. (2022). 


5.6. Carbon cycling influence 


Understanding the geodynamic evolution of the Francevillian group 
has considerable implications for unravelling processes that drove per- 
manent oxygenation of the atmosphere during the GOE (Poulton et al., 
2021). This is because the Francevillian group is a rare paleogeographic 
setting enclosing the Lomagundi Event (LE)—a time when oxygen is 
thought to have first risen permanently in the atmosphere (Kump et al., 
2011; Canfield et al., 2013; Ossa Ossa et al., 2018; Mayika et al., 2020; 
Poulton et al., 2021). Unresolved debates suggest the extreme negative 
5 Coig values of <—40 %o in the Francevillian sub-basin may reflect 
oxidative weathering of an ancient uplifted Corg reservoir deposited in 
the early stages of the LE (Kump et al., 2011; Canfield et al., 2013). 
However, no evidence exists to support an uplift event, while Weber and 
Gauthier-Lafaye (2013) invoke methanotrophic activity. 

Our data suggest that negative 5 3Corg excursion in the Francevillian 
Group was initiated somewhere at the top of the banded shales succes- 
sion, before intensifying up section in conjunction with hydrothermal 
activity and authigenic P enrichment. Considerable Corg production in 
surface waters is mirrored by increasing Corg burial in the upper 
metalliferous lithologies and negative 5 3Corg covariation (Fig. 9B, C). 
Similar observations for the oxidation of seawater derived Corg in the 
Lastoursville sub-basin (Yoshida et al., 2023), further argues against the 
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uplift-oxidative weathering of an ancient Corg model (Kump et al., 
2011). Instead, our data suggest that the upwelling of abundant dis- 
solved nutrient fluids stimulated primary productivity in the continental 
margin surface waters. The released Corg invigorated anaerobic micro- 
bial activity, with biological methane cycling and anaerobic Corg 
oxidation contributing to negative 5 Cog enrichment in the sediments, 
sulfate reduction to pyrite formation, dissimilatory Fe and Mn oxide 
reduction to Fe-Mn carbonate deposition, and fixed N reduction to 
positive 5'°N enrichment. 


6. Conclusions 


The results of this study indicate that the Paleoproterozoic France- 
villian sub-basin was hydrologically restricted and even cut off from the 
global ocean to sustain stable nutrient enrichment and exchange over 
extended periods. In these conditions, up section dissolved P enrichment 
and the proliferation of benthic microbial mat feeding grounds would 
have provided a natural nutrient-rich laboratory for biological experi- 
mentation with allopatric speciation to achieve more complex macro- 
biological proportions. This evolutionary experiment is suggested to 
have been triggered by a combination of tectonic processes that drove 
nutrient enrichment and recycling, culminating in local seawater 
enrichment with P. These observations make it possible that the 
appearance of macrofossils in Franceville may mark a unique window 
into our understanding of conditions that enabled and restricted the 
evolution and disappearance of Earth’s earliest macrobiological life 
forms. Particularly, the broadly anoxic and nutrient-poor Paleoproter- 
ozoic Ocean water mass, together with potential confinement of the 
Francevillian sub-basin, would have presented difficult conditions for 
the Francevillian biota to take a permanent and global foothold (e.g., 
Reinhard et al., 2016). When coupled to the absence of similar mani- 
festations in the mostly anoxic nutrient-deficient interval bracketing the 
Francevillian and Ediacaran biotas, our observations imply that complex 
macrobiological forms may have evolved at least twice in Earth histo- 
ry—i.e., shortly after the GOE and the NOE. These conclusions draw 
support from recent Zn systematics (El Albani et al., 2023; Ossa Ossa 
et al., 2023) and genomic dating (Craig et al., 2023) swinging the needle 
for when eukaryotic life potentially first emerged from ~1700 Ma to- 
wards ~2100 Ma. 
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